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ABSTRACT 

We have studied ~ 2100 early-type galaxies in the SDSS DR3 which have been detected by the 
GALEX Medium Imaging Survey (MIS), in the redshift range < z < 0.11. Combining GALEX 
UV photometry with corollary optical data from the SDSS, we find that, at a 95 percent confidence 
level, at least ~ 30 percent of galaxies in this sample have UV to optical colours consistent with some 
recent star formation within the last Gyr. In particular, galaxies with a NUV — r colour less than 
5.5 are very likely to have experienced such recent star formation, taking into account the possibility 
of a contribution to NUV flux from the UV upturn phenomenon. We find quantitative agreement 
between the observations and the predictions of a semi-analytical ACDM hierarchical merger model 
and deduce that early-type galaxies in the redshift range < z < 0.11 have ~ 1 to 3 percent of their 
stellar mass in stars less than 1 Gyr old. The average age of this recently formed population is ~ 300 
to 500 Myrs. We also find that 'monolithically' evolving galaxies, where recent star formation can be 
driven solely by recycled gas from stellar mass loss, cannot exhibit the blue colours (NUV — r< 5.5) 
seen in a significant fraction (~ 30 percent) of our observed sample. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
galaxies: fundamental parameters 



1. INTRODUCTION 

One of the most important unresolved debates in con- 
temporary astrophysics concerns the formation mech- 
anism of early-type galaxie s . In t he classical 'mono- 
lithic ' model (e.g. iLarsonl H97l [Ch iosi and Carrarol 
|2002| ). early-type stellar populations form in a sin- 
gle, short, highly efficient burst of star formation at 
high redshift (z >> 1), which is followed by passive 
ageing to present day. Such a scenario can explain 
many of the observed properties of early-type galax- 
ies, b oth at local and high redshift (jChiosi and Carrarol 
2002), witho ut invoking hierarchi c al me r ger-driven pro- 
cesses (e.g. iToomre and Toomrd Il972t iToomrei 11977 



Kauffmann et al.l 



19931: Somerville and Primac 



Benson et al. 1 120001: ICole et al. H2000HHatton et allhooa 
Khochfar and Burkertl 120031) . Observational evidence 
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that might conclusively rule out either the monolithic 
or the merger -based scenario has remained elusive (see 
iPeeblesI (|2002j) for a critical review of the two paradigms) . 
The properties of optical colour-magnitude relations 
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(CMRs) of early-type galaxies have often been taken 
as evidence in favour of the monolithic scenario. 
The lack of redshift evoluti on of the slope an d scat- 
ter in op t ical C MRs (e.g. iBower et all [1991 119981: 
Ellis et all H997t iStanford et all 119981: iGladders et al.1 



19981 : Ivan D okkum et al.ll2000l ) is consistent with a high 
formation redshift (z > 2), as is the evolution in their 
zero-point. However, the predicted star formation histo- 
ries (SFHs) of early-type galaxies in the merger paradigm 
make such a conclusion far less clear cut. For example, 
the predicted SFHs of (cluster) early-types in the merger 
scenario are quasi-monolithic, with an overwhelming ma- 
j ority of the stellar mass forming before a redshift of 1 
(|Kavirai et al. 2005). A careful treatment in the merger 
framework suggests that, the observed optical CMR 
can also be reconciled comfortably with merger models. 
From a photometric point of view, optical colours are, 
at best, degenerate with respect to the competing theo- 
ries of early-type galaxy formation and it is very difficult 
to discriminate betw een them using optical colours alone 
(jKavirai et al.ll2005h . 

Evidence does exist for morphological evolution in 
galaxies, suggesting that formation mechanisms of 
early-types are at least not uniquely monolithic. Al- 
though approximately 80 percent of galaxies in the 
cores of present day clusters have early- type morpholo- 
gies (jDresslerl 119801 ) , a higher fraction of spiral galax- 
ies have been reported in clusters at 0.3 < z < 

0.8 (e.g. iButcher and Oemlerlll984; iDressler et al.l ll997l: 
iCouch et all 119981 : Ivan D okkum et al J 120001) , along with 
increased rates of merger and interaction e vents (e.g. 
ICouch et all Il998t Ivan Dokkum et all 119991) . This is 
supported by recent results which suggest that the 
mass density on the red sequence (which is domi- 
nated by early-type systems) has doubled since z = 
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1 (jBell et all l2004f ). The strengths of age-sensitive 
spectral indices observed in some early-type systems 
require luminosity-weighted ages which are consistent 
with the presence of at least so me recent star forma- 
tion (RSF) in thes e s ystems (e.g. ICaldwell et al.1 120031 : 
iTrager et all I2000al|b1 : IProctor and Sansoml 120021 ) . Fur- 
ther more , somfijiejajb^^ 

5128 (e.g. IReikuba et al.l l 2001l: IPeng et al.ll2003l : lYi~et~aTl 
120041: IKavirai et al.ll2005l; IReikuba , et al.ll2004ft . NGC205 
( Hod gelll973HBurstein et allll988f). N GC 510 2 dPritchetj 
Il979t Burstem et alj 119881: IDeharveng et alj|1997l ) and 

certain early-type sys tems recently mapped by the SDSS 
(|Fukugita et al.ll2~004f ) exhibit unambiguous signs of sig- 
nificant ongoing or recent star formatio n. Others, 
such as NGC 2865 (iBica .and Alloinl Il98l. NGC 5128 
(ISchiminovich et al.Hl994ft . NGC 3921 (|Schweizer et all 
|1996|) provide direct residual evidence of the recent in- 
teractions that created them (see also papers by the 
SAURON collaboration e.g. Falcon-Barroso et al. 2005 
and references therein). It is worth noting that, con- 
trary to the traditional notion of early-type galaxies be- 
ing (cold) gas-poor systems, observational evidence over 
the last thirty years (e.g. from IRAS) have shown that 
this is not the case. Cold gas, which supplies the fuel for 
star formation, has been detected in significant numbers 
of nearby early-type systems (e.g. IKnapp et alj 119891 : 
iKnapp and Rupen ! 1996; Young 2005], also see the review 



by IKnapp! (|1999f )). Thus it is perhaps correct to expect 
RSF in early-type galaxies, rather than consider them 
inert, passively evolving systems! 

Given the apparently-evolving early-type fraction and 
accumulating spectro-photometric evidence for RSF in 
nearby early-type galaxies, it is vital to establish whether 
RSF may only be a sporadic feature of a few early-types 
or a more widespread phenomenon in the early-type pop- 
ulation. The confirmation of low- level star formation in a 
significantly large sample of low-redshift early-type galax- 
ies would put powerful constraints on models for their 
formation. 

Spectroscopic indicators of RSF are already available, 
such as the commonly used Hp index, higher order 
Balmer lines such as H 1 and Hs and the D4000 break. 
From a photometric point of view, rest-frame ultraviolet 
flux is highly sensitive to RSF. Given the apparent de- 
generacy in optical colours, ultra-violet photometry of- 
fers the best route to finding tell-tale signatures of RSF 
in early-type systems at low redshift - if indeed such star 
formation is present in these systems! 

In a recent work, iFerreras and Silkl (|2000al . hereafter 
FS2000) studied a sample of early-type galaxies in the 
cluster Abell 851 at z = 0.41. They used F300W (rest- 
frame near- ultraviolet (NUV), 2000 angstroms) and op- 
tical photometry to perform one of the first studies of 
the A?"/7 V-optical CMR in a sample of early-type galax- 
ies. They found that the slope and scatter of the NUV- 
optical CMR was consistent with some early-types hav- 
ing ~ 10 percent of their stellar mass in stars younger 
than ~ 500 Myrs. Detailed modelling of this data 
(|Ferreras et al.|[2002l ) indicated that secondary bursts of 
star formation at low redshifts lead to a natural expla- 
nation of the larg e scatter in the NUV -op tical CMR. 

More recently, IDeharveng et al.l ([2002) studied the 
rest-frame flux around 2000 angstroms from a sample of 
82 nearby early-types using data from the FOCA, SCAP 



and FAUST experiments. The focus of this study was 
primarily to investigate the far-ultraviolet (FUV) emis- 
sion of early- type galaxies, thought to originate from old 
extreme horizon tal branch (EHB) stars and their progeny 
(|Yi et al.lll997h . However, their analysis suggested that 
the 2000 — V colours in some early-type galaxies are sig- 
nificantly bluer than would be expected from old popu- 
lations alone and cannot be explained without invoking 
some RSF. 

Low redshift photometry in the far-ultraviolet (FUV; 
1530 angstroms) and near-ultraviolet (NUV; 2310 
angstroms) passbands from the GALEX mission (Mar- 
tin et al. 2005), unprecedented both in terms of its 
quality and quantity, provides a unique opportunity to 
study the RSF-sensitive UV emission from a variety of 
nearby early-type galaxies across a range of luminosities 
and environments. In this paper we study ~ 2100 early- 
type galaxies detected by GALEX, in the redshift range 
< z < 0.11. We focus mainly on the NUV passband, 
because the FUV is sensitive not only to RSF but also 
to the UV upturn flux from EHB stars which could be 
expected in old early- type populations. The NUV is less 
sensitive to UV upturn and therefore a better RSF indi- 
cator. 

A preliminary study of the UV emission of early-types, 
based on GALEX detections of SDSS early-type galaxies 
listed in the catalog of Ber nardi et al. (20 03d, B2003 
hereafter) was presented in lYi et alj (|2005h . Compar- 
ing the photometry of early-types in this sample to the 
spectral energy distribution (SED) of a strong nearby 
UV- upturn galaxy (NGC 4552), they did not find more 
than 2 galaxies (out of 162) which fit the typical UV- 
optical shape of a system with a significant amount of 
UV upturn flux. They concluded that the large scat- 
ter in the NUV — r colours cannot be generated from 
a scatter in late- stage stellar evolut ion alone. We also 
direct readers to iBoselli et al.l ((2005) who were the first 
to perform a detailed study of the UV CMR in the Virgo 
cluster. However, we note that our sample covers a wide 
variety of environments and we sample almost two mag- 
nitudes deeper (in r-b and) in dense environments (see 
Schawinski et al.ll2005f) compared to the Virgo sample of 



Boselli et a 1 .11201)5^ 



We begin this study by describing the constr uction of 
an ea rly-type catalog, similar to that used in lYi et all 
(2005), but in which the morphological classification of 
objects and removal of potentially t/U-contaminating 
AGN are performed in a more robust manner. Using a 
simple parametrisation of the SFH, we identify galaxies 
which appear to have had some star formation within the 
last Gyr. We use a semi-analytical ACDM model, cali- 
brated to accurately reproduce the (cluster ) optical CMR 
in th e redshift range < z < 1.27 (see IKavirai et al] 
120051 ). to reproduce the observed NUV CMR. Finally, we 
discuss the comparative roles of quiescent and merger- 
driven processes in driving the residual star formation 
found in some of our early-type sample and investigate 
whether an extreme monolithic scenario is capable of 
reproducing the NUV colours of large (> L*), blue 
(NUV — r < 5.5) early-type galaxies in our sample. 

2. SAMPLE SELECTION 

2.1. Constructing an early-type catalog from the SDSS 

DR3 
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B2003 were the first to construct an SDSS catalogue 
of ~ 9000 early-type galaxies, selected using a combina- 
tion of SDSS pipeline parameters (see also Bernardi et 
al. (2003b, c,d) for a comprehensive study of the prop- 
erties of galaxies in this catalogue). Their catalogue 
contains galaxies with high i-band concentration indices 
( r 5o/ r 90 > 2.5), spectral classifications typical of early- 
type galaxies, and where the deVaucouleurs fit to the sur- 
face brightness profile is more likely than the exponential 
fit. While such automated prescriptions are very efficient 
at selecting a reasonably robust early-type galaxy sam- 
ple, it suffers from two shortcomings. 

Firstly, it is apparent from visually inspecting a sample 
of Bernardi early- types that it contains late-type contam- 
inants, such as sideways-on spirals, Sa-type galaxies and 
objects where a central dominant bulge is surrounded 
by faint spiral features. While clearly not early-type 
galaxies, these objects all pass the selection criteria. The 
UV emission from early-type galaxies is predicted to be 
weak in any galaxy formation scenario. The presence of 
late-type contaminants is therefore a particular problem, 
since such galaxies will be substantially bluer than true 
early-types and our conclusions regarding the presence 
of star formation in the early-type population could be 
severely affected 10 . 

Secondly, spectral selection criteria, such as those used 
in B2003, immediately excludes any early-type galaxy 
which may contain star formatio n and therefore s pectra l 
lines, such as those reported by iFukugita et al] (2004). 
Since one of our main aims is to test galaxy formation 
models and attempt to discriminate between the mono- 
lithic and merger driven paradigms, we must take special 
care not to bias our sample towards one scenario at the 
outset! Clearly, removing early- type systems with signs 
of star formation will bias the sample towards the pas- 
sive monolithic models and make any conclusions regard- 
ing the recent star formation activity in early-types less 
meaningful. 

Since robust morphological classification inevitably re- 
quires visual confirmation of early-type candidates, our 
approach to constructing an early-type catalog relies 
more on visual inspection of each candidate than auto- 
matic extraction based on pipeline parameters. We make 
our initial selection of objects based on a single SDSS 
parameter - fracDev, which is the weight of the deVau- 
couleur's fit in the best composite (deVaucouleur's + ex- 
ponential) fit to the galaxy's image in a particular band. 
We extract all galaxies which have fracDev > 0.95 in 
g,r and i bands. The short wavelength g-band traces 
disk or spiral arm structures, while r and i trace the cen- 
tral bulge. The robustness of this procedure has been 
checked by applying it to a sample of 200 randomly se- 
lected nearby early-type galaxies from the SDSS DR3. 
We find that 90 percent of the early-types in this sample 
(confirmed by visual inspection) have fracDev > 0.95. 
This implies that, on average, this criterion would pick 
up 90 percent of the galaxies in a typical sample of SDSS 
galaxies. Reducing the value of fracDev = 0.95 would 
admit more early-types but also increase contamination, 
resulting in a larger sample of galaxies to be processed 

10 We must note, however, that given the size of the B2003 
sample, performing a visual inspection of all galaxies is clearly 
impractical! 



visually. A fracDev value of 0.95 gives the best trade- 
off between selecting a sample dominated by early-type 
objects and a sample with a significant amount of con- 
tamination. Note that we have chosen not to make a dis- 
tinction between elliptical and SO galaxies in this study 
since it is difficult to discriminate between them effec- 
tively using the SDSS images. It is worth noting, how- 
ever J _Jhatjth^sjtwo_populations could behave differently 
(e.g lBoselli end1l2005h . 

Having extracted our initial working sample of galax- 
ies, we perform a visual inspection of each early-type can- 
didate. However, our sample of galaxies spans a large 
range in luminosity and redshift. Our ability to accu- 
rately classify any galaxy depends on the clarity of its 
features on the SDSS image - this in turn depends on 
its redshift and apparent luminosity. Thus we must first 
determine the limiting redshift and luminosity at which 
we can trust our classifications. 

The average SDSS exposure time is 50 seconds. To per- 
form this check we use 24,000 second exposure COMBO- 
17 images of SDSS galaxies, at redshifts between 0.11 and 
0.13 (kindly provided by Chris Wolf). The high exposure 
time of the COMBO-17 images allows us to resolve all 
features within the galaxies in question. We find that, for 
an apparent r-band magnitude of 16.8, our classification 
based on the SDSS images matches that based on the 
COMBO-17 images. We therefore apply a redshift cut 
of 0.11 and an apparent r-band magnitude cut of 16.8 
to our sample, to extract a magnitude limited sample of 
~ 350 SDSS early-type galaxies (see ISchawinski et al.1 
(|2005f) for more details). 

2.2. Cross-matching with GALEX MIS detections 
We direc t readers to iMartin and the GALEX Team! 



(pOOBT) and iMorrissev and the GALEX Team! 1)20051 ) for 
a description of the scope and performance of the 
GALEX satellite. The magnitude-limited SDSS early- 
type sample described in the previous section is now 
cross-matched with data from 595 GALEX fields, im- 
aged in the medium depth (MIS) mode. GALEX has 
a fiducial angular resolution of 6 arcseconds. The posi- 
tional matching is performed within a more conservative 
4 arcseconds. We eliminate any GALEX objects that 
have multiple SDSS matches within 6 arcseconds because 
it is impossible to tell which companion dominates the 
UV flux detected by GALEX. We find that the fraction 
of galaxies in our magnitude limited sample (r < 16.8, 
z < 0.11) detected by GALEX is around 90 to 95 per- 
cent (Figure [TJ) . The detection rate drops sharply after 
z = 0.11. 

2.3. AGN diagnostics: optical analysis 

Since we are studying the (generally weak) rest-frame 
UV flux from early-type galaxies here, it is important 
to ascertain the level of possible contamination of the 
UV continuum due to non-thermal emission from AGN, 
which are common in early-type galaxies. iRich et al.l 
(2005), who employed a similar method to B2003 to ex- 
tract a sample of early-type galaxies from the SDSS and 
study their UV emission detected by GALEX, were the 
first to perform this additional check for AGN contami- 
nation, using emission lines measured by the SDSS. Fol- 
lowing their study, we perform a similar analysis to gauge 
the contribution from AGN to the rest-frame UV in our 
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Fig. 1. — Fraction of SDSS early-types in MIS fields which 
are detected by GALEX. 



galaxies, and eventually remove those which exhibit signs 
of hosting a strong AGN. 

Type I (i.e. unobscured) AGN are largely removed by 
using the SDSS spectral classification algorithm (by set- 
ting the QSO flag) during our initial extraction of galax- 
ies - any remaining Type Is are then removed completely 
through the optical and radio analyses we describe below. 
Type II (i.e. partly obscured) AGN can be distinguished 
from normal (star-forming) galaxies usin g the intensity 
ratio s of pairs of strong emission lines (Baldwin et al. 
I1981L BPT hereafter). Using the emission line ratios 
[OIII/H/3] and [NII/Ha] iKauffmann etal] (|2003l . K2003 
hereafter) have used a BPT type analysis to classify a 
large sample of SDSS galaxies into star-forming and Type 
II AGN (seyferts, LINERs and transition objects). We 
identify and remove potential AGN from our sample us- 
ing the criteria derived by K2003. 

Before calculating line ratios, we must correct for pos- 
sible absorption in these line regions and extract the true 
emission strengths of each line index. We first construct 
a library of 10,000 star formation histories (SFHs) in 
which each SFH is modelled by an instantaneous star- 
burst at z = 3 followed by a second instantaneous star- 
burst, which is allowed to vary in age and mass frac- 
tion. The choice of this model library is motivated by 
the fact that our galaxy sample contains only early-type 
objects and we find that such model SFHs give excel- 
lent fits to the stellar continua of galaxies in our sample. 
While studies such as K2003 have to deal with the full 
spectrum of morphological types, and thus require more 
elaborate model libraries, we tailor our model library to 
the narrow set of morphologies spanned by our sample. 

The procedure for correcting for absorption proceeds 
in the standard way, by first comparing the spectrum of 
each observed galaxy (excluding major line regions) to 
the model library and finding the best-fit model. The 
absorption strengths of the OIII, H/3, Nil and Ha lines 
in the best-fit model are then computed and subtracted 
from the measured line strengths in the observed spec- 
trum. For each galaxy which exhibits all four emission 
lines with a signal-to-noise ratio (S/N) greater than 3, 
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Fig. 2.— Plot of [OIII/H/3] vs. [NII/Ha] used to find and remove 
AGN in our sample. The curved line represents the demarcation 
between star-forming systems and Type II AGN - galaxies are con- 
sidered star-forming if they lie below this line. Red points repre- 
sent galaxies identified as early-type after visual inspection. Blue 
points are galaxies identified as contaminants in our initial sample 
of galaxies. The small grey dots are K2003's SDSS Type II AGN 
catalog. It is worth noting that the star-forming locus (see Figure 
1 in K2003 for comparison) is almost completely empty, which is 
expected since our sample is composed exclusively of early-type 
galaxies. In addition, the majority of galaxies which lie on the 
star-forming locus are contaminants. 



we construct the line ratios, place the galaxy on a plot 
of [OIII/H/3] vs. [NII/Ha] and classify it as either star- 
forming or a Type II AGN, based on the criteria derived 
by K2003. 

In FigureHwe plot [OIII/H/3] vs. [NII/Ha] used to find 
and remove AGN in our sample. The curved line rep- 
resents the demarcation between star-forming systems 
and Type II AGN - galaxies are considered star-forming 
if they lie below this line. Red points represent galax- 
ies identified as early-type after visual inspection. Blue 
points are galaxies identified as contaminants in our ini- 
tial sample of galaxies. The small grey dots are K2003's 
SDSS Type II AGN catalog. It is worth noting that the 
star-forming locus (see Figure 1 in K2003 for compari- 
son) is almost completely empty, which is expected since 
our sample is composed exclusively of early-type galax- 
ies. In addition, the majority of galaxies which do lie 
on the star-forming locus are contaminants. This optical 
analysis identifies ~ 25 percent of our early-type sample 
as Type II AGN. 

2.4. AGN diagnostics: radio analysis 

Galaxies which either do not have the four required 
emission lines, or in which emission line detections have 
S/N < 3, cannot be treated using the optical analysis 
discussed above. In most studies these galaxies are as- 
sumed to have no contribution from AGN. However, we 
opt not to make this assumption and confirm the results 
of the optical analysis for galaxies which remain unclas- 
sified using their radio emission. There are two radio 
surveys which overlap with our ea rly-type sample - the 
FIRST survey (|Becker et al.|[l995T) . which has an angu- 
lar resolution of ~ 5 arcseconds and a completeness limit 
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Fig. 3. — Comparison of radio luminosities of galaxies classified 
by the optical analysis as normal (crosses) and Type II AGN (all 
other symbols, see legend for specific AGN classifications.) Type 
II AGN tend to dominate this plot above a radio luminosity of ~ 
W 22 WHz~ 1 , while most of the normal (i.e. star-forming) objects 
lie below this threshold. 



of 1 mjy, and the NVSS survey, which has an angular 
resolution of 45 arsceconds and a completeness limit of 
2.5 mjy. 

In this study we use the FIRST survey due to its supe- 
rior resolution and depth - we find that including NVSS 
in our analysis does not add to or alter our conclusions. 
The high resolution of FIRST poses a slight problem 
for sources where the radio emission is not very cen- 
tralised. For example, galaxies with radio lobes have 
multiple, spatially resolved detections in FIRST. We 
therefore cross-match our early-type sample with FIRST 
detections using a large search radius of 30 arcseconds. 
If there are multiple detections within this search radius 
we add the fluxes of all the detections corresponding to 
the early-type object in question. 

To identify a luminosity threshold above which radio 
activity could be attributed to AGN, we compare the 
radio luminosities of galaxies classified by the previous 
optical analysis as normal (i.e. star-forming) to the lu- 
minosities of those classified as Type II AGN. We show 
this comparison in Figure[3] Crosse show normal galaxies 
while other symbols correspond to various classes of Type 
II AGN. We find that Type II AGN tend to dominate this 
plot above a radio luminosity of ~ 10 22 W H z _1 , while 
most of the normal (i.e. star-forming) objects lie below 
this threshold. We use this threshold and remove any 
galaxy which remains unclassified in the optical analysis 
but shows a radio luminosity above lO^WHz^ 1 . This 
removes ~ 3 percent of the sample of galaxies which sur- 
vive the optical analysis. 

We must note, however, that galaxies with a high ra- 
dio luminosity do not necessarily host AGN. There is a 
degeneracy between radio emission from AGN and from 
genuine star formation which is not quantified accurately 
by using a single threshold value. By applying this cut 
we might indeed exclude interesting early-types galaxies 
which are genuinely star-forming. However, since the ra- 
dio analysis only affects a small fraction of our sample, 



TABLE 1 

Summary of criteria used in constructing the 
SDSS-GALEX early-type catalog. 



Criterion Reason 

r < 16.8 Robust morphology 

(from COMBO-17 comparison) 
z < 0.11 Robust morphology, 

detection rate stable at 90 percent 
Emission line analysis AGN removal for emission line 
galaxies with S/N > 3 
L(radio) > 10 22 WHz~ 1 Further removal of possible AGN 
not classified by the emission 
line analysis 



we include it in our sample selection, as it would only 
make our sample slightly more conservative in terms of 
star-forming early-type galaxies. 

2.5. The final catalog 

Table [TJ summarises all the criteria used in the con- 
struction of our final catalog of early-type galaxies - the 
final sample contains 2116 galaxies. 

Figure[4]shows the optical (top panel) and the NUV—r 
(bottom panel) colour-magnitude relations (CMRs) of 
our early- type galaxies. Galaxies are colour-coded ac- 
cording to their redshifts. The small grey crosses rep- 
resent galaxies which are rejected based on the optical 
or radio analyses described in the previous sections. At 
first glance, the tight optical CMR typical of early-type 
populations is in stark contrast to the the large spread in 
colours evident in the NUV CMR. While the scatter in 
the optical CMR is approximately 0.05 mag, the spread 
in the NUV colours is almost 6 mags with a (bi-weight) 
scatter of approximately 1 mag. The properties of the 
NUV CMR immediately lead to two conclusions. 

Firstly, it seems inconsistent with a 'monolithic' evo- 
lution scenario for all galaxies in our sample. If all early- 
types were dustless simple stellar populations (as is fre- 
quently assumed from their tight optical CMRs in the 
context of the monolithic picture), we might expect the 
NUV CMR to be just as tight as its optical counter- 
parts. In reality, any age dispersion, coupled with the 
presence of varying amounts of UV upturn, might be 
expected to induce some scatter in the NUV CMR. A 
lower limit to the NUV — r colour in such a hypotheti- 
cal monolithic CMR can be estimated by looking at the 
NUV — r colour of the nearby strong IW-upturn galaxy, 
NGC 4552. We use a composite spectrum of NGC 4552, 
in which the UV spectral ranges are construc ted from 
IUE a nd HUT data (see caption of Figure 1 in lYi et all 
(fl99l for more details on sources in the construction of 
this composite SED), to estimate this lower limit. This 
composite SED yields a NUV — r colour of 5.4. We 
should note, however, that the UV component, derived 
from IUE and HUT data, samples the UV bright core 
of NGC 4552, so the value of 5.4 is likely to be bluer than 
the NUV — r colour averaged over the entire galaxy. This 
makes NUV — r — 5.4 an optimistic estimate for the blue 
end of a monolithic CMR. However, in a monolithic sce- 
nario all the galaxies would have NUV — r > 5.4, which 
is clearly not the case for the galaxies in our sample. 
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Fig. 4. — The optical (top panel) and the NUV — r (bottom 
panel) colour-magnitude relations (CMRs) of our early-type galax- 
ies. We intentionally show the optical CMR on the same scale as 
its UV counterpart, to highlight the significant difference in their 
respective scatters. Galaxies are colour-coded according to their 
redshifts. The small grey crosses represent galaxies which are re- 
jected based on the optical or radio analyses. 



We can therefore assume that the excess scatter in the 
NUV CMR is due, at least in part, to some recent star 
formation. 

Secondly, it is apparent that the red envelope is not 
very well-defined i.e. we do not see a sharp red edge to 
the NUV CMR as is frequently seen in optical CMRs. 
This is likely to be caused by varying levels of dust 
(and therefore gas) in these galaxies which smear out 
the tightness in the red sequence. 

Note that NUV colours are a factor of 4 more sensi- 
tive to dust than optical colours such as (g — r), so the 
presence of a small amount of dust can have a propor- 
tionately larger effect in the NUV than in the optical 
spectral ranges. 

In Figure [5] we compare the BPT rejected galaxies to 
galaxies in our final catalog. The BPT rejected galaxies 
prefer bluer NUV — r colours, which is not unexpected 
since the continuum AGN contribution makes the NUV 
colour bluer. There appears to be no bias with redshift 
which fits our expectation that the fraction of AGN host- 
ing galaxies will not change significantly over such a small 
redshift range. In Figure EJ we show the NUV luminosi- 
ties of early- types in our sample as a function of redshift. 
The yellow region shows the detectable M^uv space, 
assuming the fiducial detection limit of the MIS fields 
(m(NUV) = 23, dashed line). The absolute NUV mag- 
nitude of the detected early-types and their errors are 
shown by the blue points. Also shown is the expected 
absolute NUV magnitude of a passively evolving giant 
early-type galaxy (My — —23) which forms at z = 3, 
assuming solar (solid line) and twice solar (dashed line) 
metallicities. Since the monolithic evolution shown is for 
an extremely bright early-type galaxy, detections above 
this line immediately imply that, for at least some early- 
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Fig. 5. — Comparison of BPT rejected galaxies to galaxies in 
our final catalog. The BPT rejected galaxies prefer bluer NUV — r 
colours, which is not unexpected since the AGN contribution to 
the continuum would be expected to make the NUV colour bluer. 
There appears to be no bias with redshift - we do not expect the 
fraction of AGN hosting galaxies to change significantly over this 
redshift range. 



types, a monolithic SFH is not consistent with their de- 
tected NUV flux. 

As shown in Figure [TJ approximately 10 percent of the 
SDSS early-types in the GALEX field-of-view remain un- 
detected. To explore the properties of these undetected 
galaxies, we construct, for each non-detection, a syn- 
thetic absolute NUV magnitude, based on the assump- 
tion that the galaxy is dustless, has solar metallicity and 
is formed at z = 3. We plot these synthetic magnitudes 
as red points in Figure [51 The filled red circles corre- 
spond to undetected galaxies with r < 16.8 and the small 
red points are non-detections with r > 16.8. We find 
that, while some non-detections are close to the detec- 
tion limit, a significant number of these galaxies should 
be detected, under the assumptions used to construct 
their synthetic magnitude. However, both a super-solar 
metallicity and/or a small dust content can be sufficient 
to push these galaxies below the detection limit - we indi- 
cate this, in Figure O by showing the predicted deviation 
of the NUV luminosity from the synthetic magnitude as- 
sumptions for a small dust content, e.g. E(B — V) = 0.1, 
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Fig. 6. — SDSS early-types detected by GALEX: the absolute 
NUV magnitude has to be in the yellow region for the galaxy 
to be detected by GALEX. The detection limit corresponding to 
the MIS fields is shown by the dashed curve. The absolute NUV 
magnitude of the detected early-types and their errors are shown 
by the blue points. Also shown is the expected evolution of the 
absolute NUV magnitude of a monolithic giant early-type galaxy 
(My = —23) which formed at z=3 and has solar (solid line) and 
twice solar (dashed line) metallicity. 



and non solar metallicities. 

3. QUANTIFYING THE RSF IN THE GALEX EARLY- TYPES 

The aim of this section is to quantify the RSF in our 
early-type sample. In particular, we would like to estab- 
lish how many of our galaxies have photometry consistent 
with any star formation within the last Gyr, at the 95 
percent confidence level. 

Before we begin, it is instructive to consider the scope 
of the NUV in computing age determinations. Figure [7] 
shows the intrinsic NUV — r colour of a simple stellar 
population (SSP) as a function of age and metallicity 
(solar: dotted line, half-solar: solid line and twice-solar: 
dashed line). The evolution in the NUV — r colour slows 
considerably after 1-2 Gyrs, as the NUV flux from a 
recent starburst decays significantly in this timescale - 
the NUV is therefore an effective age indicator for only 
~ 1 Gyr after the starburst in question. In other words, 
we can expect the UV (+ optical) photometry of our 
galaxies to be able to distinguish only between galaxies 
which have had a very recent starburst, within the last 
1 Gyr and galaxies that havent. The NUV, however, is 
blind to any star formation which is much older than 1 
Gyr old. 

We fit the observed colours of the GALEX early-types 
across the UV and optical spectral ranges to colours pre- 
dicted by a library of model SFHs, in which an initial 
starburst at high redshift (z = 3) is followed by a second 
starburst which is allowed to vary in mass fraction and 
age between z = 3 and present day. Both starbursts are 
assumed to be instantaneous. Each model in the library 
has two main parameters: tyc, the age of the second 
starburst, in the range 10~ 2 to 15 Gyrs and fyc, the 
mass fraction of stars formed in the second starburst, in 
the range 10 -4 to 1 ('YC stands for Young Component). 
We explore models with metallicities Z/Zq in the range 
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Fig. 7. — NUV — r (top) and u — r (bottom) colours of a simple 
stellar population (SSP) as a function of age. The solid line shows 
SSPs with half solar metallicity, the dotted line represents solar 
metallicity and the dashed line indicates twice solar metallicity. 
The evolution in the NUV — r colour slows after 2 Gyrs because 
the NUV flux from a recent starburst decays significantly in this 
timescale. 



0.5 to 2.5 and intrinsic E(B — V) values in the range 
to 0.15. We use the standard dust prescription given 
in ICalzetti et al.l (120001) to compute the reddening for a 
given value of E(B — V). The stellar models used to 
co mpute the model library in this study (|Yi et al.lll997h 
lYll2003f ) adopt the universal Salpeter IMF. For galaxies 
without FUV detection, we use a four colour (NUV — r, 
g — r, r — i and r — z) comparison, while for those with 
FUV detection we use five colours (FUV — r, NUV — r, 
g — r, r — i and r — z). We perform a \ 2 test with the 
appropriate degrees of freedom by minimising the sum of 
the normalised residuals 
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where C™ and C° are the n model and observed 
colours respectively and a n is the uncertainty in the 
residual C™ od - C° bs . The average errors in the SDSS 
optical magnitudes are approximately 0.01 mag. The av- 
erage uncertainty in the GALEX FUV and NUV mag- 
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nitudes are ~ 0.25 mag and ~ 0.15 mag respectively. 
The uncertainties in stellar models are taken to be ap- 
proximately 0.05 mag for optical pa ssbands and 0.1 mag 
for the UV passbands (see lYilr2003[ ). Since, our observed 
sample contains galaxies at various redshifts, models are 
redshifted to z = for each comparison, thus avoiding 
uncertainties due to K-corrections. The derived model 
fits are therefore in the rest-frame of each observed early- 
type galaxy. 

We use this procedure to derive best-fit values of 
(tYCifyc) an d % 2 likelihoods in (tycfyc) space for 
each galaxy. Our entire sample can be broadly divided 
into three categories - we illustrate this by showing the % 2 
contour maps of four example galaxies in Figure [5J The 
X 2 minimum is marked with x and we indicate the 68 per- 
cent and 95 percent confidence contours. Galaxies 1 and 
2 (top row) have clear signatures of RSF, since their 95 
percent contours are contained completely within tyc < 
1 Gyr. Galaxy 3 shows a large age-mass degeneracy 
which appears in most cases because the NUV flux from 
a small young starburst can be indistinguishable from 
that due to a larger but older starburst. Degenerate 
contours can also be the result of large errors in the ob- 
served photometry which allow many models to fit the 
data with values of (reduced) \ 2 < 2. Therefore, even 
though Galaxy 3 has a best-fit tyc around 1 Gyr, we 
cannot rule out the possibility of younger or older pop- 
ulations being present, purely from its UV and optical 
photometry. Finally, Galaxy 4 is likely to have old or 
intermediate age populations, since its 95 percent confi- 
dence contour lies exclusively at ages above 2 Gyrs. 

We note that, although they used a similar model for 
their analysis, the age-mass degeneracies in fits to our 
sample of galaxies are smaller than those studied by 
FS2000. This is partly due to the fact that optical un- 
certainties in SDSS measurements are ten times smaller 
than those used by FS2000. Together with the additional 
FUV constraint (FS2000 only used rest-frame NUV), 
this allows a much smaller part of the (fyc>tYc) space 
to fit the colours of the observed early-types. 

In Figure [9] we summarise the best-fit values of 
(tyc, fyc) for our entire sample of galaxies. Red points 
indicate the best-fit values for early-type galaxies in our 
sample while orange points indicate best-fit positions of 
contaminants (galaxies which are rejected by the visual 
inspection or rejected by the AGN diagnostics). We find 
an increased scatter of contaminants towards lower val- 
ues of tyc (i.e. younger ages). We also show the 95 per- 
cent confidence contours of a galaxy with strong, moder- 
ate and weak UV upturn (UVX) flux. We generate the 
strong UVX contour using the observed SED of NGC 
4552. The moderate and weak UVX contours are com- 
puted by reducing the UV flux (< 3500 angstroms) in the 
NGC 4552 SED by a factor of 2 and 6 respectively. In 
Figure fTUl we show the distribution of best-fit tyc values 
in our sample of galaxies - the inset shows the distribu- 
tion of values restricted to the range < < 1 
Gyrs. Figures [H and [TU] suggest that the UV + opti- 
cal photometric properties of the majority of early-type 
galaxies in our sample are not compatible with a single 
rapid burst at high redshift followed by passive evolu- 
tion. Such a conclusion is difficu lt to achieve with optical 
CMRs studied in the past (e.g. iBower et al.lll992L 119981 : 
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Fig. 8. — Four examples of y contour maps. The \ minimum 
is marked with 'x' and we indicate the 68 percent (lo - ) and 95 
percent confidence contours. Galaxies 1 and 2 (top row) have clear 
signatures of RSF, since the 95 percent contours are contained 
completely within tyc < lGyr. Galaxy 3 shows a large age-mass 
degeneracy because the NUV flux from a small young starburst 
can be Galaxy 4 is clearly likely to have old or intermediate age 
populations since its best-fit tyc ls greater than 3 Gyrs and its 95 
percent confidence contour does not extend far below 2 Gyrs. 



Ellis et alJ 119971: IStanford et alj 1 19981 iGladders et~ai] 
1998t Ivan D okkum et al.ll2000f ) . which are all consistent 



with purely passive monolithic scenarios. 

We now classify our galaxies using their \ 2 contours, 
as follows: if the 95 percent confidence contour lies ex- 
clusively below a tyc of 1 Gyr then we label the galaxy 
as young (e.g. galaxies 1 and 2 in Figured]). If the 95 
percent contour spans ages above and below 1 Gyrs, the 
galaxy is classified as degenerate (e.g. galaxy 3 in Fig- 
ure |HJ • Finally, if the 95 percent contour lies exclusively 
above an age of 1 Gyr then the galaxy is classified as 
old (e.g. galaxy 4 in Figure [8]). As indicated in Figure 
[9j some galaxies have best-fit values which fall within 
the UVX contours defined by the observed SED of NGC 
4552. These galaxies clearly have SEDs which resemble 
that of typical UV upturn galaxies but they may be ini- 
tially labelled as young. However, since it is impossible 
to determine whether the shape of their SED is a result 
of RSF or UVX or a combination of both, we re-label 
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Fig. 9. — Best-fit values of fyc and tyc- 
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Fig. 10. — Distribution of best-fit tyc values in our sample of 
galaxies - the inset shows the distribution of values restricted to 
the range < tyQ~ %t < 1 Gyrs 



these galaxies as degenerate. 

In Figure [Til we plot the u — r and NUV — r CMRs 
for the GALEX early-types and colour-code each galaxy 
using this classification scheme. We find that the optical 
CMR is unable to distinguish between different levels of 
RSF, since the various classifications cannot be distin- 
guished in the optical CM space. However, the NUV 
CMR shows a clear layered structure. The galaxies clas- 
sified as young (blue symbols) are bluest in the NUV — r 
colours and clearly separated from the old locus (red 
symbols) and degenerate class of galaxies (green sym- 
bols) forming an intermediate layer. Blue galaxies reside 
almost exclusively below NUV — r ~ 5.5 (for comparison, 
the NUV -r colour of NGC 4552 is 5.4). Degenerate and 
old galaxies lie almost exclusively above NUV — r ~ 5.5. 

We note that this method indicates the presence of 
RSF in a very general sense. Two instantaneous star- 
bursts are clearly an oversimplification of the true SFH 
of early-type galaxies. The position of a galaxy in 
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Fig. 11.— The optical u — r CMR (top panel) and the NUV 
CMR (bottom panel) for the GALEX early-types colour-coded us- 
ing the following classification scheme: if the 95 percent confidence 
contour does not extend to ages higher than 2 Gyrs then the galaxy 
is classified as young (e.g. galaxies 1 and 2 in Figure |8j; if the 95 
percent confidence contour spans ages below and above 2 Gyrs the 
galaxy is classified as degenerate (e.g. galaxy 3 in Figure[8j; if 95 
percent confidence contour lies exclusively above an age of 2 Gyrs 
then the galaxy is classified as old (e.g. galaxy 4 in Figure [5J. 
Young galaxies are shown in blue, intermediate galaxies are shown 
in green and old galaxies are shown in red. Note that the photome- 
try has not been K-corrected because this requires the assumption 
of a template - it is evident that an old passively evolving popu- 
lation is not good template for the local early-type population, at 
least for UV photometry. 



(fYCjtyc) space is a function of the shape of its SED 
and our efforts in this section have centred on approxi- 
mating the shape of the true SED of each galaxy (using a 
library of fyc and tyc) through colour fitting. Although 
the quantities fyc and tyc do not have a direct physical 
meaning, the two-starburst method probes whether there 
has been any star formation in these systems within the 
last Gyr. Since our classifications take account of the 95 
percent confidence contours, i.e. the part of parameter 
which enclose 95 percent of the probability, we are able 
to derive conclusions about the presence of such RSF at 
the 95 percent confidence level. Having identified the 
signature of RSF in an object, we must then appeal to 
the various formation paradigms (see Sections 5 and 6 
below) to motivate the star formation that generates the 
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observed RSF. 

4. THE BLUE FRACTION 

Galaxies classified as young and shown in blue in Fig- 
ure [IT] are systems with a clear signature of RSF. In 
Figure [T^] we show a stacked likelihood map of all blue 
galaxies, constructed by adding the likelihood maps of in- 
dividual blue galaxies (c.f. top row of Figure[8]) together. 
The black region indicates the part of parameter space 
where the average likelihood is higher than 95 percent. 
The dark and light grey regions correspond to regions 
where the probability is between 70 and 95 percent and 
30 and 70 percent respectively. Thus, the black region 
represents the values of fyc and tyc which are, on aver- 
age, most likely for the blue population. The most likely 
value of fyc is ~ 1 percent and tyc is ~ 0.5 Gyrs. 

In Figure [13] we split the NUV CMR into narrow red- 
shift bins and indicate the blue population in each bin. 
It is apparent from Figure [13] that, as a result of con- 
structing a (apparent) magnitude limited sample, we are 
not sampling identical parts of the luminosity function 
at different redshifts. Comparison of the blue fraction 
across our redshift range is therefore meaningful only for 
the part of the luminosity function which is sampled com- 
pletely at all redshifts - looking at the highest redshift 
bin this appears to be for galaxies with M(r) < —21.3. 

In Figure [14] we plot the evolution of the blue frac- 
tion with redshift. Recall that ~ 10 percent of early- 
type galaxies in the SDSS sample remain undetected by 
GALEX (Figure [6]) . The true value of the blue fraction 
should take these galaxies into account. We therefore 
assume that none of the non-detections are blue and ad- 
just the blue fraction accordingly. This is a fairly strong 
assumption because the undctectability of these galaxies 
is a function of both their dust content and metallicity. 
It is possible that some of these galaxies do have star 
formation but are heavily dust-reddened and/or metal- 
rich. Nevertheless, without their UV photometry this is 
a question we cannot answer. 

The solid curve in Figure [14] takes into account all 
galaxies in each redshift bin, while the dotted curve con- 
siders only galaxies with M(r) < —21.3. Note also that 
at very low redshifts (z < 0.05), SDSS has a bright detec- 
tion limit - spectroscopy is not possible for bright galaxies 
(M(r) < —22.7) which are extremely close (z < 0.05). 
Hence we do not properly sample the luminosity func- 
tion in the first two redshift bins, although the blue 
fraction should remain unchanged due to this undersam- 
pling. The red curve in Figure [T4l indicates that the frac- 
tion of blue systems for bright galaxies (M(r) < —21.3) 
varies between ~ 27 and 34 percent with an average value 
of ~ 30 percent. Within the errors there is negligible 
evolution in the blue fraction within the redshift range 
< z < 0.11. Nevertheless, both these points make our 
blue fraction a conservative estimate. 

It is worth noting that the blue fraction computed in 
this study is a lower limit on the fraction of early-type 
galaxies which contain RSF. It is possible that some of 
the galaxies classified as degenerate also contain a certain 
level of RSF, but this is impossible to verify given their 
photometry and/or observational uncertainties. In ad- 
dition, we have assumed that undetected ellipticals are 
not blue although it is clearly possible that they have 
some star formation but are heavily dust-reddened or 
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Fig. 12. — Stacked likelihood map of blue galaxies, constructed 
by adding the likelihood maps of each blue galaxy together. The 
black region indicates the region where the average likelihood is 
higher than 95 percent. The dark and light grey regions correspond 
to regions where the probability is between 70 and 95 percent and 
30 and 70 percent respectively. The black region represents the 
values of fyc an d tyc which are, on average, most likely for the 
blue population. 



very metal-rich. This is, however, unquantifiable and 
any correction to the blue fraction could be expected to 
be negligible within errors. 

5. COMPARISON TO THE MERGER PARADIGM 

Since we have established that RSF does exist in a sig- 
nificant fraction of the early-type population, it is natu- 
ral to ask whether the properties of our local early-type 
sample can be reproduced in the currently popular hi- 
erarchical merger paradigm. We therefore explore the 
predicted NUV CMR in the semi-analytical framework 
using the GALICS model, which combines large scale N- 
body simulations with analytical recipes for the dynam- 
ical evolution of baryons within dark matter haloes. We 
direct readers to lHatton et al.l (2003, H2 003 hereafter) 
for sp ecifics regarding the fiducial model. iKavirai et al.1 
( 2005) have used this model to accurately reproduce the 
the observed optical CMR in cluster early-type galaxies 
in the redshift range < z < 1.27 (see their Section 
2). We use the blind UV predictions from this optically- 
calibrated model to study the colours of the GALEX 
early-type sample. Note that galaxy morphologies in the 
model are defined by the ratio of S-band luminosities 
of the disc and bul ge components, which correlates well 
with Hubble type (jSimien and de Vaucouleurj|1986f ). A 
morphology index is defined as 

7 = exp-L B /L D (2) 

so that pure discs have 7 = an d pure bulges have 
1=1. Following lBaugh et al.1 (|1996f ). early- type galaxies 
have / < 0.507 11 . 

We approach the problem from the point of view of 
reproducing the (NUV — r) colour distribution using an 
ensemble of model galaxies in the redshift range spanned 

11 In this scheme ellipticals have / < 0.219, SOs have 0.219 < 
/ < 0.507 and spirals have / > 0.507. 
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Fig. 13.— The evolution of the NUV CMR with redshift. Galax- 
ies classified as young in the colour-fitting analysis are shown using 
filled circles. 
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Fig. 15.— Evolution of the UV flux with age. TOP: Relative 
UV fluxes of simple stellar populations. BOTTOM: Relative mag- 
nitudes of simple stellar populations shown in the top panel. The 
dotted lines indicate the peak throughput of the FUV and NUV 
filters. All models assume solar metallicities and are taken from 
the public BC2003 distribution. 
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Fig. 14. — Evolution of the blue fraction with redshift. The 
solid curve takes into accounts all galaxies in each redshift bin (see 
Figure 1130 , while the dotted curve considers only galaxies with 
M(r) < —21.3. Note that, due to arguments given in Section 4, 
the blue fraction presented here is a conservative estimate. 



by the observed sample. In traditional (optical) CMR 
studies, the colour distribution is often parametrised by 
the slope and scatter of the correlation. This is useful 
because the optical correlations are typically tight i.e. 
have small scatter. But since the NUV CMR has more 
structure than can be expressed through the slope and 
scatter alone, matching only the slope and scatter within 
errors (which are large due to the significant dispersion 
in colours) is no longer sufficient. 

5.1. Treatment of very young stars 

Since we are modelling the UV spectrum in this study, 
we must address some additional issues regarding the flux 
from (very) young stars in the model. The evolution of 
UV flux with the age of the J71/-emitting young stars 
is highly non-linear. Figure [15] shows the evolution of 
the UV spectrum with age, on both flux and magnitude 
scales for a solar meta l licity stellar model from the public 
iBruzual and CharlotJ (|2003l ) distribution. It is apparent 
that the UV flux of a 1 Myr population is ~ 4-5 orders 
of magnitude greater than that of a 1 Gyr population. 
Like other models of galaxy formation, star formation in 
GALICS is essentially proportional to the cold gas avail- 
able in the galaxy - stars are therefore produced continu- 
ously. In addition, the GALICS model provides high age 
resolution (1 Myr) for very young stars. Therefore in a 



12 



S. Kaviraj et al. 



significant number of galaxies, a residual fraction of very 
young stars (< 10 Myrs old) is present. 

There are two main issues regarding such young stars. 
Very young stars are typically embedded in birth clouds 
(BCs). During this embedded phase the attenuation due 
to dust in the ambient ISM is compounded by the dust 
present in the outer HI envelopes of the BCs. Thus for 
embedded stars, the combined optical depth due to the 
birth cloud (BC) and the inter-stellar medium (ISM) can 
be several times large r than due to the ISM alone (e.g. 
iCharlot and Fallll2000D . While most studies suggest that 
birth clouds are tr ansient features with lifetimes of a few 
tens of Myrs (e.g. Bli tz and Shul 11980b iHartmann et al.l 
1200 1|) , a recent work (|Scoville and Wirsonll2004h suggests 
that the low efficiency of star formation in these clouds 
implies significantly longer lifetimes (~ 100 Myrs). Re- 
cently formed stars may still migrate early from such 
long-lived BCs, since the BCs feel hydro-dynamic forces 
while the stars do not. However, the low gas density in 
early-type galaxies imply lower hydro-dynamical forces, 
which could mean that migration timescales are longer 
(Nick Scoville, private communication). 

In this study we apply a differential extinction scheme, 
in which the extinction applied to flux from embedded 
stars is a factor /x greater than the extinction felt by stars 
which are not embedded - calculated from the mean H 
column density (see H2003, Section 6.2). The key in- 
gredients in this scheme are the factor /j, and the age 
limit T max of stars which are considered to be embedded. 
Following the studies mentioned above, we experiment 
with various values of T max in the range 10-100 Myrs 
and a range of values of ix b etween 3 and 10. We note 
that ICharlot and Falll (|200Q[ ) suggest < fi >~ 3, albeit 
with a large scatter at lower and higher values - how- 
ever this is mainly for starbursting galaxies which have 
significantly different star formation histori es than the 
early -type galaxies in our sample (see also iKong et alJ 
l200l . 

A second issue regards the modelling of star formation 
itself. Star formation is modelled as a continuous pro- 
cess proportional to the mass of cold gas in the system 
- as long as there is gas to fuel star formation, there are 
always some newborn i.e. zero-age stars predicted in the 
model at any epoch. However, the question is whether, 
when galaxies are observed, true zero-age stars or indeed 
very young stars are actually seen. In reality, star forma- 
tion is unlikely to be an idealised continuous process but 
probably occurs in random bursts of varying strengths 
and durations. Galaxies are therefore observed sometime 
in between two starbursts i.e. there is a small time delay 
between a galaxy's most recent starburst and the point at 
which it is observed. For example, NGC 205, an 'actively ' 
star-forming dwarf elliptical (e.g. Burstein et alJ |1988|: 
Wilcots et all 119901: iDorman et all Il995t iBertola et all 
1995ft in which the two most rec ent starbursts occur red 
~ 20 Myrs and ~ 500 Myrs ago (|Wilcots et al.lll990D . is 
noticeably b luer in UV-optical co lours than other ellipti- 
cal galaxies (|Burstein et al.|[l988f h This suggests that it 
may not be correct to model all early-type galaxies with 
very young stars (< 20 Myrs old). 

If we assume that starbursts in real galaxies are ran- 
domly, i.e. Poisson distributed, in time, then the dis- 
tribution of time delays between starbursts is exponen- 



tial 12 . The shape of the exponential is determined by 
the number of starbursts per Gyr (k). In addition, the 
time lag distribution must be truncated at some value 
Imax, since the time lag should generally be small and 
could have a minimum value l m in (e.g. l m in = 20 Myrs 
and Imax = 500 Myrs, from the NGC 205 arguments 
above). In this study we explore the (Kd m i n d max ) pa- 
rameter space for values of l m in in the range 0-0.05 Gyrs 
and Imax in the range 0.2-0.5 Gyrs. We assume that 
k = 10 i.e. that a starburst has, on average, a timescale 
of ~ 0.1 Gyrs, so that galaxies tend to have ~ 10 bursts 
per Gyr. 

5.2. Reproducing the NUV CMR in the merger 
paradigm 

A one-to-one comparison between the model and ob- 
served photometry requires that the two samples are con- 
sistent in terms of the magnitude and redshift ranges 
covered by each. We therefore take the following points 
into account while performing our comparisons: 

• From Figure [14] we see that, due to the magnitude 
limit of our observed sample, luminosity function 
coverage is complete only for galaxies brighter than 
M(r) = —21.3. We therefore restrict our com- 
parison only to model galaxies with M(r) < —21.3. 



• As mentioned before, the SDSS spectroscopic 
sample has a bright limit - bright galaxies which 
are very close do not have spectroscopy and 
therefore do not appear in our observed sample. 
A large sample of galaxies drawn from the SDSS 
DR3 indicates that this bright limit lies around 
M(r) — —22.7 for z < 0.05. We therefore remove 
such galaxies from the model sample. 



• The GALEX detection limit (m(NUV) < 23) 
implies that some galaxies predicted in the model 
will not be 'GALEX-detectable'. The comparison 
must therefore include only model galaxies which 
are GALEX-detectable. 



• A fourth point which does not particularly affect 
our comparison but is worth noting is that the 
model itself has a completeness limit (H2003, 
Section 8.2 and Table 3), because small haloes 
are not fully resolved in the simulation due to the 
mass resolution of GALICS 13 . The completeness 
limit for GALICS in r-band is around -20.4. 
While the effects of such a lack of resolution 
may propagate over the formal completeness 
limit of -20.4, this does not affect our compar- 
ison if we restrict ourselves to M(r) < —21.3, 

12 This is a general property of Poisson distributions 

13 In GALICS the merging history is driven by the N-body 
backbone on which galaxies are 'painted'. The mass resolu- 
tion is therefore determined by the mass of the dark matter 
particles (~ 10 10 Mp,} . Other models, (e.g. IBenson et al.l 120001 ; 
IKhochfar and Burkert] l2003h . use the Prcss-Schcchtcr formalism 
(Press and Schcchtcr 1974), with theoretically infinite resolution, 
so that the merging history of even the smallest objects is fully 
resolved 
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FlG. 16. — Comparison of observed and model photometry in 
the case where fi = 10 and T max = 50 Myrs, i.e. the BC + ISM 
extinction is ten times that of ISM alone and BCs survive for 50 
Myrs. The open red circles represent the observed photometry. 
All other colours represent model galaxies at various redshifts (see 
legend in plot). 



since this is well within the completeness threshold. 



• Finally, all model SEDs are redshifted to present 
day before the comparison, to avoid uncertainties 
due to K-corrections. This is important, because 
without robust low redshift templates for the UV 
spectra of early-type galaxies, we cannot reliably 
estimate UV K-corrections - we have already 
established that a purely passively evolving tem- 
plate does not match at least 30 percent of the 
early-type galaxies at low redshifts. Later in this 
section we present UV K-corrections based on the 
photometry of our best-fit models. 



We generally find that reproducing the observed pho- 
tometry without invoking timclags is difficult. The as- 
sumption of moderate BC survival spans (~10-50 Myrs) 
leaves the model photometry with roughly correct scat- 
ter but with a zero-point which is consistently too blue, 
regardless of how heavily attenuated the flux from young 
stars is. To illustrate this point we show, in Figure [TBI 
the case where fi — 10 and r max — 50 Myrs, i.e. the BC 
+ ISM extinction is ten times that of the ISM alone and 
BCs survive for 50 Myrs. Note that, for clarity, in this 
and other similar figures that follow, we draw a random 
sample of model galaxies which has the same size as the 
observed sample. While the top panel of Figure [16] indi- 
cates that the scatter is reasonably reproduced the colour 
distributions in the bottom panels clearly show that the 
zero-point of model photometry is clearly too blue by ~ 
0.5 mag in such a scenario. 

The situation can be significantly improved by invok- 
ing longer BC lifetimes close to 100 Myrs - Figure [TH] 
shows the comparison between the model and observed 
results when an exponential BC lifetime distribution, 
highly skewed towards an upper limit of 100 Myrs is 
used, along with a BC extinction of 5 times that due 
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FlG. 17. — Exponential distribution of BC lifetimes employed to 
produce the comparison shown in Figure [18] 




WW 




1.0 
o.s 

0.6 
0,4 
0,2 
0,0 



. ■ Model 






■ Observed 








i 






- 


') 2 + 


6 


5 



FlG. 18. — Comparison of observed and model photometry in the 
case where an exponential distribution of BC lifetimes is employed 
(see Figure [TTt . The extinction assumed due to the BCs is 5 times 
that due to the ISM alone i.e. fj, = 5. The open red circles repre- 
sent the observed photometry. All other colours represent model 
galaxies at various redshifts (see legend in plot). 



to the ISM alone. We find that this comparison gives a 
reasonable fit to the photometry, although as mentioned 
before, it is unclear whether recently formed stars within 
BCs linger within them for the lifetime of the cloud itself. 

Invoking small timclags into the model photometry, 
while keeping BC lifetimes short, produces similar re- 
sults. The observed photometry can be reproduced best 
if we assume a minimum timelag of 20 Myrs (c.f. NGC 
205 arguments given above), a maximum timelag of ~ 
200-300 Myrs and an exponential timelag distribution 
constructed on the assumption that galaxies have, on 
average, 10 starbursts per Gyr. We also assume that 
stars less than 30 Myrs old reside in BCs and experi- 
ence three_Jhnes_the extinction due to the ISM alone 
(c.f. iCharlot and F all (2000). Using the terminology of 
the previous section this translates to: /j = 3, T max = 30 
Myrs, k = 10, l mm = 20 Myrs and l max = 200 Myrs. 



14 



S. Kaviraj et al. 




HW-r NUV-r 

Fig. 19. — Comparison of observed and model photometry in the 
best-fi t model where fi = 3, r max = 30 Myrs (c.f. Chariot and Fall 
(EES), l min = 20 Myrs and l max = 200 Myrs. The open red circles 
represent the observed photometry. All other colours represent 
model galaxies at various redshifts (see legend in plot). 
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Fig. 20.— Evolution of the JV(7 V CMR in the best-fit model. The 
red dots indicate galaxies which are below the detection threshold 
of GALEX. The grey dots lie in the region below the complete- 
ness limit of the GALICS model. The detection limit encroaches 
severely on the model sample at the upper end of our redshift range 
= 0.11). 



Note that the BC + ISM extinction only affects stars 
which are between 20 and 30 Myrs old since, by con- 
struction, there are no stars younger than 20 Myrs in 
this model (since lmin = 20 Myrs). This comparison is 
shown in Figure [IT)] Henceforth, we refer to this sce- 
nario as the 'best-fit model'. As might be expected, the 
goodness of fit is more sensitive to the minimum timelag 
lmin - scenarios where l m in < 20 Myrs make both the 
zero-point too blue and the scatter too large to fit the 
observations. The fit is less sensitive to the maximum 
timelag but we use l ma x = 200 Myrs, the lowest value of 
lmax which gives a good fit. 

As the bottom panels of Figure [TH] indicate, the best- 
fit scenario reproduces the colour distribution of the ob- 
served sample well. The evolution of the (complete) 
NUV CMR in this scenario is shown in Figure |2"0"1 Note 
that the number of galaxies in each redshift bin is around 
~ 6000. The red dots indicate galaxies which are below 
the detection limit of GALEX i.e. they are not GALEX- 
detectable. The grey dots lie in the region below the 
completeness limit of the GALICS model. The detection 
limit encroaches severely on the model sample at the up- 
per end of our redshift range (z = 0.11). 

While it is apparent that the best-fit model reproduces 
the UV photometry of the observed sample, we now ap- 
ply the two-component analysis used in Section 3 to the 
predicted photometry from the best-fit model. The posi- 
tion of an object in {fyc, Wc) space is a direct function 
of the shape of its SED, which in turn drives its photo- 
metric properties both in the UV and the optical spectral 
ranges. A comparison, in this parameter space, between 
the positions of the model galaxies and the observed pop- 
ulation indicates the overall robustness of the fit, not only 
in terms of the NUV colour but in terms of the shape 
of the complete SED. The derivation of (fyc, Wc) f° r 
model galaxies uses the same process as that used for the 
observed sample. We show this comparison in Figure [2"T1 
Note that we show all available galaxies in the observed 
sample and similarly all model galaxies above the GAL- 



TABLE 2 

The average properties of recent star 
formation in model early-types. 



z 


</> 


< T m > 


<T V > 


< T r > 


0.03 


0.014 


0.61 


0.55 


0.55 


0.07 


0.021 


0.62 


0.55 


0.57 


0.11 


0.025 


0.50 


0.37 


0.38 



Note. — < / > is the average 'RSF' 
mass fraction i.e. the stellar mass frac- 
tion within the last Gyr. < t x > is the 
average age of this RSF fraction, in Gyrs, 
weighted by the quantity indicated in the 
subscript x. Subscript m indicates mass, 
v indicates F-band luminosity and r in- 
dicates r-band luminosity. 



ICS completeness limit - we do not restrict the analysis 
to galaxies with M(r) < —21.3. We find that there are 
no systematic differences between the positions of model 
and observed galaxies in this parameter space, which im- 
plies that the observed photometry across the UV and 
optical spectral ranges is reproduced reasonably well by 
the best- fit model. Table [2] summarises the properties of 
the recent star formation in the galaxies in the best-fit 
model. On average, low-redshift early-type galaxies are 
predicted to have ~ 1.5-2.5 percent of their stellar mass 
formed within the last Gyr. The mass-weighted age of 
this star formation is between 0.4 and 0.6 Gyrs and its 
luminosity-weighted average age is around 0.5 Gyrs. 

We end this section by computing K-corrections 
(Hog g et al.ll2002l ) from our best-fit models. As suggested 
by our previous analysis, purely old, passively evolving 
populations are clearly not adequately good templates 
for the UV spectra of early-type galaxies, even though 
the UV emission in these galaxies is weak. In Figure [221 
we show K-corrections for UV colours from our best-fit 
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Fig. 21. — The two-component method applied to model early- 
types in the GALICS hierarchical merger model which is used in 
this study. Black points represent model galaxies and the red points 
represent the observations. Note that this comparison is not re- 
stricted to galaxies with M(r) < —21.3. All observed galaxies and 
all model galaxies above the completeness limit of M(r) = —20.4 
are plotted. There are no significant systematic differences be- 
tween the positions of model and observed galaxies in this param- 
eter space, implying that the observed photometry across the UV 
and optical spectral ranges is reproduced reasonably well by the 
models. 
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model in the redshift range spanned by our observed sam- 
ple (black dots with error bars). We also compare these 
K-corrections to those obtained from a passively evolving 
9 Gyrs old SSP (blue) and the elliptical t emplate from 
the Kinney spectral atlas of galaxies (red, iKinnev et alJ 
119961) . The use of a 9 Gyr SSP is motivated by the fact 
that Bernardi et al. (2003d) find that their early-type 
sample can be best-fit by a 9 Gyr passively evolving SSP, 
formed at high redshift. 

5.3. The role of mergers and progenitor bias 

Star formation in the merger paradigm proceeds 
through one of two channels - quiescent star formation 
governed by the mass of available cold gas and the dy- 
namical timescale of the galaxy (H2003, Section 4.1) or 
merger driven star formation as a result of dynamical 
interactions between galaxies (H2003, Section 5). As 
Figure [23] indicates, it is possible for galaxies to have 
experienced very recent mergers, within the last Gyr of 
look-back time. The top panel in Figure [53] shows the 
redshifts of the most recent mergers experienced by early- 
type galaxies at z = 0. These are essentially the dynam- 
ical ages of these galaxies i.e. the epoch at which the 
galaxy in its present (morphological) form was created. 
The bottom panel shows the average trend in dynamical 
ages as a function of luminosity (size) . In a hierarchical 
picture of galaxy formation larger early-types are assem- 
bled later, i.e. have smaller dynamical ages, although the 
ages of their stellar populations are, on average, older 
([Kaviraj et alJl2005f ). 

We investigate the impact such recent mergers may 
have on the colour of the early-type remnants, and the 
extent to which they may affect the blue colours exhib- 
ited by some early-type galaxies. In Figure [2"4l we com- 
pare the UV colour distribution of galaxies which have 



Redshift 

Fig. 22. — K-corrections for UV colours from our best-fit model 
in the redshift range spanned by our observed sample (black) . Also 
shown are K-corrections obtained from a passively evolving 9 Gyrs 
old stellar model (blue) and the elliptical template from the Kinney 
spectral atlas of galaxies (red). 



experienced recent mergers (within the last Gyr), with 
the global colour distribution of all early- type galaxies. 
We find that recent mergers do not have a preponderance 
of blue galaxies, suggesting that mergers at low redshift 
are predominantly dry and not the cause for the blue scat- 
ter in the UV colours of early-type galaxies - star forma- 
tion from the quiescent mode is at least as important as 
that due to mergers. 

We finish this section by briefly mentioning the is- 
sue of progenitor bias, which affects all studies of early- 
type galaxies at z > 0, in the framework of the merger 
paradigm. Fundamental to merger-driven morphological 
transformations is the implication that, at higher red- 
shifts, a progressively larger fraction of mass that even- 
tually resides in present-day e arly-types is potentially 
locke d up in late-type systems ([van Dokkum and Franxl 
2001). Like all other early- type studies, this investiga- 
tion has attempted to trace the evolution of the local 
early-type population by looking only at their early-type 
progenitors in the redshift range < z < 0.11 - any non 
early-type progenitors have been excluded, for the simple 
reason that it is not possible to easily identify them! Al- 
though we expect progenitor bias to be smallest at low 
redshifts (indeed this is the big advantage of studying 
nearby early-type galaxies!), it is instructive to quantify 
this effect and gauge its possible impact on the blue frac- 



16 



S. Kaviraj et al. 




0,001 1 i r, ,.*, ,' 

-24 -23 -22 -21 -20 -19 -IB 



14 




n l , , , , i , , , , i , , , , i , , , , i , , , , i , , , , i , , , , 
-23.0 -22.5 -22.0 -21.5 -21.0 -20.5 -20.0 -19.5 
<M(V)> 

Fig. 23. — TOP: Redshifts of most recent mergers experienced 
by early-type galaxies at z = 0. These are essentially the dynam- 
ical ages of these galaxies i.e. the epoch at which the galaxy in 
its present (morphological) form was created. BOTTOM: Average 
dynamical ages of galaxies as a function of their luminosity (size). 
As expected, in a hierarchical picture of galaxy formation, larger 
galaxies are assembled later i.e. have smaller dynamical ages, al- 
though the ages of t heir stellar populations are, on average, older 
flKavirai et alj|2003) . 



tion derived in Section 4. 

The top panel in Figure 1251 indicates the extent of pro- 
genitor bias by plotting the fraction of the z = carly- 
type galaxies which have already undergone their last 
merger and assumed their final early-type morphology 
as a function of redshift. At z ~ 0.1, roughly 10 percent 
of present-day early-type galaxies are still 'in pieces'. If, 
for example, these 'pieces' were blue late type systems, 
then we could have underestimated the blue fraction in 
Section 4. However, we have already shown (Figure [24]) 
that recent mergers in this redshift range have a simi- 
lar UV colour distribution to the rest of the population 
i.e. that recent mergers are not preferentially bluer. In- 
deed, mergers at low redshift have a high probability 
of involving at least one early-type progenitor (bottom 
pane l of Figure [2"5"|) . Note that our results are consistent 
with iKhochfar and Burkertl (|2003f ) who were the first to 
present a detailed analysis of the progenitor morpholo- 
gies of early- type galaxies. Thus, since progenitor bias is 
small within our redshift range and mergers tend to be 
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Fig. 24. — LEFT: Comparison between the UV colours of early- 
type galaxies which have had recent mergers (blue) within the last 
Gyr of look-back time with the colour distribution of all galaxies. 
RIGHT: The UV colours of early-type galaxies as a function of 
its dynamical (last-merger) age. The dy nam ical age is explained 
in the text and in the caption to Figure 1231 We find that recent 
mergers do not have a preponderance of blue galaxies, suggesting 
that mergers at low redshift are predominantly dry and not the 
cause for the blue scatter in the UV colours of early-type galaxies. 



dry, neglecting non-early-type progenitors will not affect 
the blue fraction derived in this study. 

6. REVISITING THE MONOLITHIC HYPOTHESIS 

The monolithic hypothesis remains the simplest and 
perhaps most elegant explanation for a remarkably 
wide range of pro perties of large early-type galaxies 
(e.g. iPeeblesI 12002) . Indeed, a careful treatment us- 
ing N-body/SPH simulat ions (jChiosi and Carraroll2002t 
iTantalo and Chiosil I2002D suggests that monolithic for- 
mation can reproduce a wide spectrum of chemo- 
photometric properties of large early- types. Until now 
these properties have been restricted to the optical spec- 
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Fig. 25. — TOP: The fraction of local (z = 0) early-type galaxies 
which have already attained their final (early-type) morphology, 
i.e. had their last merger, as a function of redshift. BOTTOM: 
Progenitor types in binary mergers. At low redshift mergers have 
a h igher probability of having a t least one early-type progenitor. 
See Khochfar and Burkcrt (2003) for a detailed analysis of the pro- 
genitor morphologies of early-type galaxies. 



trum - the UV photometry presented in this paper pro- 
vides, for the first time, a new constraint on the forma- 
tion of such large early-type systems. While it is appar- 
ent that a merger-based treatment is able to reproduce 
the UV + optical photometry of local early-types fairly 
accurately, the RSF mass fractions (~ 1 to 3 percent) 
involved are small. The presence of young stars in these 
galaxies indicates that there is fuel, i.e. cold gas, avail- 
able for star formation at low redshift. It is, therefore, 
natural to ask whether the required level of fuel could be 
supplied in a traditional monolithic scheme. 

Monolithic models propose a high star formation effi- 
ciency at high redshift, which drives a violent episode of 
star formation resulting in the bulk of the galaxy form- 
ing within a short timescale at high redshift (e.g. z < 2). 
However, stars created in the primordial burst would re- 
cycle a fraction of their mass back into the ISM through 
stellar winds and supernova ejecta. A fraction or all of 
this internally recycled gas could fuel further star for- 
mation, without invoking any external gas accretion or 
interactions. The aim of this section is to investigate 



whether, given reasonable assumptions for the evolution 
and chemical enrichment of such a monolithic system, 
we can reproduce the RSF and the observed UV pho- 
tometry of large early-type systems that form part of the 
blue fraction 14 . The reason for restricting ourselves to 
large blue galaxies is because, as mentioned before, red 
galaxies can be reconciled with a monolithic scenario, es- 
pecially if we assume that part of the scatter in the red 
sequence is determined by varying levels of intrinsic dust 
in these galaxies. This monolithic analysis is useful be- 
cause, if low level star formation seen in the observed blue 
population can be produced simply by returned material 
from stars formed in the primordial burst, then there 
would, in principle, be no need to invoke a merger based 
scenario to produce early- type galaxies! 

6.1. Simulating 'monolithic collapse' 

The model used here (see iFerreras and Silkl [2000b) 
follo ws the standa rd chemical enrichment equations 
/e.g. iTinslevi 1 198ft. Stellar yields are taken from 
iThielemann et alj (l!996f) for stellar masses M* > l OM^ 
and from Ivan den Hoek and Groenewegenl (|1997f ) for 
lower mass stars. We use a standard Salpeter IMF and 
note that assume instantaneous recycling but not instan- 
taneous mass loss, so that the lifetime of stars of various 
masses is taken into account. 

We simulate a scenario where (gaussian) gas infall, 
with a short timescale (~ 1 Gyr), occurs at 'high' red- 
shift. Most studies which support the monolithic hypoth- 
esis agree on a formation epoch somewhere before z — 2 
- it is difficult to pinpoint it any further using the avail- 
able data. Our conclusions are not sensitive to the exact 
redshift of formation so we use a fiducial value of z — 3. 
The 'formation redshift' is taken as the epoch at which 
the gas infall peaks. The star formation rate is govern ed 
by the generally accepted Schmidt law (|SchmidtJll959l ): 

* = Cg a (3) 

where ^ is the star formation rate, C is the star forma- 
tion efficiency, g is the gas mass and a = 1.5. 

We explore the behaviour of the system using two free 
parameters. The first and most significant free parame- 
ter is the star formation efficiency C in Equation [3] Low 
values of C result in more extended star formation, while 
high values of C result in faster consumption of gas and 
a shorter episode of star formation - note that C is not 
numerically restricted to a value between and 1. As 
expected, increasing C to arbitrarily high values does 
not result in the star formation increasing without limit, 
simply because the system cannot consume gas faster 
than it is being put in! We explore the behaviour of the 
system in this 'limiting' high C regime, where star for- 
mation essentially tracks the gas infall very closely i.e. as 
soon as gas becomes available, it is almost instantly con- 
verted to stars. As a result, the overwhelming bulk of the 
galaxy is already constructed by the time gas infall ends 
at high redshift - mimicking the fundamental feature of a 
'monolithic collapse'. However, star formation does not 
stop when gas infall ceases. Gaseous ejecta returned from 

14 While we could simply explore the amount of recycled mate- 
rial expected in a typical monolithic simulation, it is necessary to 
simulate the accompanying chemical enrichment because the NUV- 
r colour is a function of both the total amount and the metallicity 
of the recycled gas 
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Fig. 26. — LEFT-HAND COLUMN: Star formation rate (top row), gas infall (middle row) and the evolution of star and gas (bottom 
row) for B fixed at and various values of C (see legend in the top row). RIGHT-HAND COLUMN: Average quantities - rows from top 
to bottom show average age, metallicity, alpha-enhancement, fraction of stars formed within the last Gyr (the 'RSF' fraction), u — r colour 
and NUV — r colour respectivel y, as a function of C (with B fixed at 0). Also indicated, in the relevant rows, are the average values of 
[m/H] and [Mg/Fe] computed by Tragcr ct al. (2000a) for local early-type galaxies. 



dying stars is gradually released back into the ISM over 
a Hubble time, providing potential fuel for further star 
formation. It is the quantity of this recycled gas that we 
investigate here - in particular, we explore whether there 
is enough gas in this recycled material to fuel late stage 
star formation and reproduce the RSF mass fractions 
seen in large blue early-types in our observed sample. 

The second free parameter is a simple prescription for 
galactic winds, which is parametrised by assuming that 
a certain fraction (B) of gaseous ejecta is permanently 
lost from the system. B varies between and 1. We are 
primarily interested in scenarios where B is small, since 
we are exploring large early-type galaxies which should 
have potential wells deep enough to retain most of the 
ejected gas. B = corresponds to essentially closed-box 
evolution after the initial gas infall stage. 

6.2. RSF in 'monolithic' early-type galaxies 

We begin by looking at the simplest scenario, where 
B = 0, i.e. all ejecta are retained by the system - this 
is also, of course, the scenario which maximises the late 
stage returned gas fraction. The B = scenario is sum- 
marised in Figure (551 The left-hand column shows the 
star formation rate (top row), gas infall (middle row) and 



the evolution of stars and gas (bottom row) for various 
values of C (see legend in the top row), with B fixed 
at 0. We note that the 'limiting' high C regime, which 
mimics monolithic collapse, is achieved for values of C 
above ~ 100. Above this value, increasing C does not 
affect the SFR (top row), because the SFR is effectively 
limited by the rate of gas infall, since stars cannot be 
produced faster than gas is being deposited in the sys- 
tem. As mentioned before, increasing C allows the SFR 
to track the gas infall function more closely, as can be 
seen in Figure [26] - the yellow SFR curve, which repre- 
sents a system with high C, peaks almost coevally with 
the gas infall, whereas the black SFR curve, which rep- 
resents the lower C regime, produces more extended star 
formation and has its peak significantly displaced from 
the peak of the gas infall. Our focus in this section will 
be exclusively on the high C regime (C > 100), in which 
star formation after the initial gas infall stage is fuelled 
only by recycled gas. 

The right-hand column in Figure [26] shows average 
quantities produced by such a scenario. The rows from 
top to bottom show the average age, average metallic- 
ity, average alpha-enhancement, fraction of stars formed 
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Fig. 27.— Same as Figure l26l except for B = 0.3. LEFT-HAND COLUMN: Star formation rate (top row), gas infall (middle row) and 
the evolution of star and gas (bottom row) for B fixed at and various values of C (see legend in the top row). RIGHT-HAND COLUMN: 
Average quantities - rows from top to bottom show average age, metallicity, alpha-enhancement, fraction of stars formed within the last 
Gyr (the 'RSF' fraction), u — r colour and NUV — r colour respectively, as a function of C (with B fixed at 0.3). Also indicated, in the 
relevant rows, are the average values of [m/H] and [Mg/Fe] computed by Tragcr ct al. (2000a) for local early-type galaxies. 



within the last Gyr (the 'RSF' fraction), u~r colour and 
NUV — r colour respectively, as a function of C (with 
B fixed at 0). We also indicate, in the relevant rows, 
the average values (and t heir formal errors) o f [m/H] 
and [Mg/Fe] computed by iTrager et all (|2000aD for lo- 
cal early- type galaxies 15 . Figure [26l indicates that, in 
the B = scenario, 'monolithic' values of C (C > 100) 
predict present-day RSF fractions of < 0.005. In com- 
parison, we find that large (> L*), blue (NUV — r < 5.5) 
galaxies in the merger models have present-day RSF frac- 
tions of ~ 0.03. The monolithic RSF fraction, driven 
purely by the recycled material from stars, is generally a 
factor 5-6 lower than what is required, in a merger frame- 
work, to reproduce the colours of large, blue early-type 
galaxies in the observed sample. 

It is instructive to explore the full range of values for B, 
although our exclusive interest in large galaxies implies 
that values of B should be low. As we show below, very 
high values of B result in excessive metal loss produc- 
ing galaxies which are too metal-poor to fit the observed 

15 Note that ITrager et al.l d2000al1 find milde r values of (fm/Hl) 
and ([Mg/Fe]) than in previous studies (e.g. [Weiss ct al] 119951 : 
ITrager et aDll9^7t IGreegidlT997h 



metallicities of local early-types. To illustrate the trend 
in B we show, in Figure 1271 the case where B = 0.3 i.e. 
30 percent of stellar ejecta are permanently lost from the 
galaxy. We find, in this case, that galaxies are gener- 
ally more me t al-poor than the mean relations derived by 
ITrager et all (|2000a|) - increasing B simply exacerbates 
this situation. 

We summarise our investigation of monolithic models 
in Figure |2"51 which plots the NUV — r colour (top left), 
RSF fraction (bottom left), mean metallicity (top right) 
and mean [Mg/Fe] (bottom right) as a function of B 
and C. We find that, in the high C regime (C > 100), 
which we take to be the traditional definition of mono- 
lithic collapse (see arguments above), models which fit 
the metallicity and alpha-enhancement expected in lo- 
cal early-type galaxies generally do not produce the blue 
colours of large blue early-types. Models with B < 0.3 
reproduce the full spectrum of average metallicities and 
alpha-enhancements observed in large early-type galax- 
ies. However, the NUV — r colour predicted by these 
models is typically above ~ 6, because RSF fractions 
contributed by stars within the last Gyr (which dominate 
the UV flux) are too low (< 0.005), compared with blue 
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Fig. 28. — Summary of average quantities as a function of B and C: NUV — r colour (top left), RSF fraction (bottom left), mean 
metallicity (top right) and mean alpha-enhancement (bottom right). 




Age (Gyrs) 

Fig. 29. — Comparison between a 'monolithic' SFH and the 
average SFH of large (> L*), blue (NUV — r < 5.5) early-type 
galaxies predicted by the merger model (which fit the colours of 
corresponding observed early- type galaxies). The top panel shows 
the non-cumulative SFH while the bottom panels compares the 
cumulative SFH. 



galaxies predicted in merger models (~ 0.03). To illus- 
trate this, we compare, in Figure [29l a 'monolithic' SFH 



(B = 0, C = 100) to the average SFH of large (> L*), 
blue (NUV—r < 5.5) early-type galaxies predicted in the 
merger model. The top panel shows the non-cumulative 
SFH, while the bottom panel compares the cumulative 
SFH. Due to the inadequate RSF fractions, blue colours 
can only be produced in a traditional monolithic scheme 
by invoking high values of B which leads to high metal 
loss and a metal-poor galaxy. However, such metal-poor 
models do not fit the observed metallicities of large early- 
type galaxies in the nearby universe. 

It is, of course, conceivable that the expelled gas, 
parametrised by B, does not escape the galaxy poten- 
tial well completely but is re-accreted by the galaxy over 
the dynamical timescale of the galaxy. However, recall- 
ing that stellar mass which dominates the rest-frame UV 
flux is created within the last Gyr of lookback time, we 
find that for reasonable dynamical timescales, expelled 
gas falling back into the galaxy will have a negligible 
effect. We illustrate this using a simple analytical ar- 
gument. Assuming that the entire expelled fraction (B) 
falls back into the galaxy and that this 'reservoir' of ini- 
tially expelled gas empties 'exponentially' we have 

ff (t)=Bexp(-*/r) (4) 
where g(t) is the expelled gas contained in the reservoir 
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and r is the dynamical timescale of the galaxy. There- 
fore, between times t — ti and t = t% the amount of gas 
falling into the galaxy is given by 

g(At) = S[exp(-ti/r) - exp(-t 2 /r)] (5) 

To make a significant impact on the UV colour of the 
galaxy, this infalling material must contribute ~ 3 per- 
cent of the stellar mass of the galaxy. Assuming that 
the infall of expelled gas starts promptly after the main 
mass of the galaxy is in place, i.e. ~ 5 Gyrs after the be- 
ginning of star formation (see Figure [26]) . the gas which 
produces the UV emitting stars falls in between t = 7 
and t = 8; assuming the Universe is ~ 13 Gyrs old. 
Hence, UV emitting stars less than 1 Gyr old will form 
a stellar mass fraction approximately equal to 

S[exp(-7/r) - exp(-8/r)] (6) 

which is required to be ~ 0.03, to produce the blue 
observed colours. However, it is difficult to achieve such 
mass fractions - reasonable dynamical timescales close to 
what is expected for large early-type galaxies (~ 1 Gyr) 
provide stellar fractions which are two orders of magni- 
tude lower than the required ~ 3 percent. The stellar 
fraction is maximum (~ 1.5 percent) for r ~ 7 Gyrs 
which is unrealistic for large early- type galaxies. Robust 
observational constraints on the amounts of recycled gas 
will aid the type of analysis presented in this section. 
Although recent studies (e.g. lYoun"gll2005h . which com- 
pare the distribution of specific angular momenta of gas 
and stars, do indicate that at least some of the gas in 
early- type galaxies is created from stellar mass loss, it 
is still unclear what fraction of the detected gas may 
have an external origin. In particular, the gas fractions 
must eventually be correlated with the intergrated UV 
fluxes from these systems, to determine the contribution 
of internally sourced gas to the UV luminosity of the 
early-type galaxy in question - such a study is currently 
in progess using GALEX and SAURON data and results 
may become available in the near future (Martin Bureau, 
private communication) . 

However, based on the analysis presented in this sec- 
tion we conclude that monolithic evolution, where RSF 
is driven solely by recycled gas from stellar mass loss, 
is not a viable channel for the production of large blue 
early- type galaxies. 

7. CONCLUSIONS 

We have studied ~ 2100 early-type galaxies in the 
SDSS DR3 which have been detected by the GALEX 
medium depth (MIS) survey, in the redshift range < 
z < 0.11. The early- type sample has been selected 
through careful morphological inspection, with poten- 
tially [/^-contaminating AGN removed through the use 
of both optical spectral data (from the SDSS) and ra- 
dio data (from the VLA FIRST survey). At a 95 per- 
cent confidence level, at least ~ 30 percent of early-type 
galaxies in this sample have optical and UV photometry 
consistent with some recent star formation within the 
last Gyr. 



Our analysis indicates that, while optical CMRs can- 
not distinguish between early-type galaxies which have 
had star formation within the last Gyr and those which 
havent, the NUV CMR is an excellent diagnostic of RSF 
- any galaxy with NUV ~r < 5.5 has a high likelihood of 
containing recent star formation (RSF), even after tak- 
ing into account the possibility of a contribution to the 
NUV spectrum from UV upturn flux. 

Comparison of the observations to predictions of a 
semi-analytical ACDM hierarchical merger model yields 
good quantitative agreement (across the UV and optical 
spectrum), if we assume that (a) very young stars (< 30 
Myrs old) are contained in birth clouds which increase 
dust extinction by a factor of ~ 3 compared to the ISM 
alone and (b) that star formation is driven by random 
(pseudo-instantaneous) starbursts which are Poisson dis- 
tributed in time, so that there is a small timelag between 
the last starburst and the observation of the galaxy - 
the timelags are distributed exponentially in time, with 
maximum and minimum values of 200 and 20 Myrs re- 
spectively. 

Combining our parametric analysis of the UV + op- 
tical photometry of the observed early-types, with the 
properties of the predicted population in the semi- 
analytical model, we conclude that early-type galaxies 
in the redshift range < z < 0.11 are likely to have ~ 
1 to 3 percent of their stellar mass in stars less than 1 
Gyr old. The (U-band luminosity weighted age) of this 
recent star formation is ~ 300 to 500 Myrs. 

Finally, we find that monolithically evolving galaxies, 
where RSF can be produced solely from recycled gas 
due to stellar mass loss, do not exhibit the blue colours 
(NUV— r < 5.5) seen in some large early-type galaxies in 
our observed sample. While the degeneracy between the 
monolithic and merger paradigms cannot be broken for 
red early-types even with UV photometry, a monolithic 
scenario is a very unlikely channel for the evolution of 
large blue early-type systems. Such blue galaxies require 
additional fuel for star formation which must necessarily 
have an external origin. 
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